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I I Abstract. We describe a first attempt to derive properties of the regular and turbulent Galactic magnetic field 

^ . from multi-frequency polarimetric observations of the diffuse Galactic synchrotron background. A single-cell-size 

r-«>. ' model of the thin Galactic disk is constructed which includes random and regular magnetic fields and thermal and 

l/^ , relativistic electrons. The disk is irradiated from behind with a uniform partially polarized background. Radiation 

1/^ ■ from the background and from the thin disk is Faraday rotated and depolarized while propagating through the 

^^ ' medium. The model parameters are estimated from a comparison with 350 MHz observations in two regions at 

'>^ , intermediate latitudes done with the Westerbork Synthesis Radio Telescope. We obtain good consistency between 

JN ■ the estimates for the random and regular magnetic field strengths and typical scales of structure in the two 

^^ I regions. The regular magnetic field strength found is a few /iG, and the ratio of random to regular magnetic field 

r^ . strength Bran/Breg is 0.7 ± 0.5, for a typical scale of the random component of 15 ± 10 pc. Furthermore, the 

CL(' regular magnetic field is directed almost perpendicular to the line of sight. This modeling is a potentially powerful 

JL I method to estimate the structure of the Galactic magnetic field, especially when more polarimetric observations 

^ ■ of the diffuse synchrotron background at intermediate latitudes become available. 
-4— » ■ 

j^ , Key words. Magnetic fields - Polarization - Techniques: polarimetric - ISM: magnetic fields - ISM: structure - 

• • ' Radio continuum: ISM 
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Vh 1- Introduction netic field along the line of sight Bf. i?M [rad m 



a 



0.8lJne[cm-^]B«[nG]ds[pc]. 

Since the first interpretation of small-scale structure in the 

,. , 1-j J- i: -uv. r^ T -u- 1.4- Use of the diffuse Galactic synchrotron emission as 

linearly polarized component oi the Galactic synchrotron -^ 

11 ju-jiuj i. i.- nuT- • 4-8- tracer of the Galactic magnetic field is complemen- 

background as being due to l<araday rotation (Wiermgaet n i i i r 

1 1nrlo^ 1, 4.- r 4.1, /^ 1 4-- i, 4- tary to the magnetic field estimates derived from pul- 

al. 1993), many observations oi the Galactic synchrotron ■' , „ „ a^, o r.i , ^ 

.., 1 •4.-J.J.4. • 1-4-- sars (e.g. Rand & Kulkarni 1989, Ohno & Shibata 1993) 

emission have shown intricate structure m polarization on ^ '^ ' ' 

1 r, • J u 4- 4- • 4- 4- 1 and extragalactic radio sources (e.g. Simard-Normandin & 

many scales, often unaccompanied by structure m total ° ^ ° 

A 14-1, 1, -4- 1, 1, • J 4-1, 4- 4-1. T? J Kronberg 1980, Clegg et al. 1992, Brown et al. 2003), in 

power. Although it has been recognized that the Faraday & j as > ;i 

,,. J, T ■ J- r .-i i-j 14- that the diffuse background can provide a continuous field 

rotation and depolarization ot the polarized synchrotron ^ ^ 

emission are due to small-scale fluctuations in the Galactic "^ ^^^ °^ ^^^^^^ ^^^^ *^<= ^^^^ ^^^^ '^^^^ *° *^^ resolution 

magnetic field, thermal electron density and/or the line ^^ ^^^ observation. Therefore this is a unique method to 

of sight, a quantttattve description of this structure has i^f^"" ^"^^^^^ ^^^ amplitudes of fluctuations in the Galactic 

proven to be extremely difficult. This is because depo- ° ■^' 

larization arises both along the line of sight and in the The RM of the synchrotron background has been 

plane of the sky (i.e. within the telescope beam), whereas used to determine the nature of distinct objects (Gray 

the rotation measure RM is an integral over the line of et al. 1998, Uyaniker & Landecker 2002, Haverkorn et 

sight of thermal electron density He and Galactic mag- al. 2003a), to estimate the uniform component of the 

Galactic magnetic field (Haverkorn et al. 2003b) or the 

Send offprint requests to: M. Haverkorn magnetic field strength and structure in supernova rem- 
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Auriga 


Horologium 


(l,b) 


(161°, 16°) 


(137°, 7°) 


size 


7°x9° 


7°x7° 


FWHM 


5.0'x6.3' 


5.0'x5.5' 


pointings 


5x7 


5x5 


bandwidth 


5 MHz 


5 MHz 


frequencies 


341, 349, 355, 


360, 375 MHz 



Table 1. WSRT polarization observations in the constel- 
lations Auriga and Horologium. Given are position and 
size of each region, the resolution, the number of point- 
ings used to mosaic the region, the frequency bandwidth 
and the central frequency in each band. 



nants (Uyamker et al. 2002), and to infer the polarization 
horizon (Uyamker et al. 2003). However, to the authors' 
knowledge, this is the first attempt to derive the turbu- 
lent component of the magnetic field (not associated with 
any discrete structure) from the diffuse synchrotron back- 
ground. 

In this paper, we present a simple model of the Galactic 
thin disk as a synchrotron emitting and Faraday-rotating 
medium, consisting of cells with a certain electron den- 
sity and magnetic field. We compare the model predic- 
tions with observed properties of the linear polarization in 
two fields observed with the Westerbork Synthesis Radio 
Telescope (WSRT), to derive estimates for several phys- 
ical parameters of the warm ISM. In Sect. |21we discuss 
the WSRT polarization observations that will be used to 
compare to the model. Section |3| describes depth depolar- 
ization in a layer that contains both synchrotron-emitting 
and Faraday-rotating material. We discuss in Sect. ^ the 
ingredients of a model for a thin disk, combined with a 
thick disk or halo providing a constant polarized back- 
ground. In Sect. we will describe the model in some 
detail, and how observational constraints can be used to 
derive estimates for parameters like the magnetic field 
strength and direction. In Sect. El we apply the model to 
our observations and discuss its results. Finally, we present 
a summary and conclusions in Sect.d 

2. The observations 

Two fields of observation in the constellations Horologium 
and Auriga, described in detail in Haverkorn et al. (2003a, 
2003b) are used to estimate parameters of the depolar- 
ization model. Observations out of the Galactic plane are 
used to avoid discrete objects like supernova remnants and 
Hii regions, which would skew the statistical information 
of the radiation that we use. Relevant properties of the 
observations are given in Tabled Undetected large-scale 
components in Stokes Q and U are not thought to be im- 
portant in these fields around 350 MHz (Haverkorn et al. 
2004). 

Fig. n] shows the linearly polarized intensity P = 
y'Q^+lP of the two regions in grey scale. The structure 
in P is uncorrelated with total intensity /, which does not 
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Fig. 2. Histograms of (from top to bottom) Q, U and / for 
5 frequencies, and RM for Auriga (left) and Horologium 
(right). Data of Q, U and / are 5 times oversamplcd, and 
only reliably determined RMs are included. In the solid 
line histogram of RM in the Auriga region, the RM gra- 
dient over the region is subtracted; the dashed line gives 
the histogram of the observed RM including the gradient. 



show any structure on scales visible to the interferometer 
down to noise level. Rotation measures RM were derived 
from the polarization angle = 0o + RMX^, where the 
ambiguity in (j) over n 180° has been taken into account 
(Haverkorn et al. 2003b). RM maps are shown as circles 
in Fig. n The RM in the Auriga region (left) shows a 
gradient of about 1 rad m^^ per degree in the direction of 
position angle 6 — —20° (N through E). 

The linear (/)(A^)-relation can be destroyed by depo- 
larization, which yields incorrect RM values. Therefore, 
we only consider "reliably determined" RAI values, where 
"reliable" is defined by (a) the reduced x^ of the linear 
0(A^)-relation Xred < 2, and (b) the polarized intensity 
averaged over frequency (P) > 20 mJy/beam (~ 4 - 5a). 

Histograms of the distributions of Stokes parameters 
Q, U and /, and oi RM are given in Fig.[21for Auriga (left) 
and Horologium (right). In the / map, point sources were 
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Fig. 1. RM maps of the observed regions in Auriga (left) and in Horologium (right), overlaid on polarized intensity at 
349 MHz in grey scale. White denotes high polarized intensity, with a maximum of 95 mJy/beam. RMs are denoted 
by white circles, and filled circles are positive RMs. The diameter of the symbol represents the magnitude of RM, 
where the scaling is given in rad m^^. Only RMs for which (P) > 5cr and reduced x^ of the linear 0(A^)-relation < 2 
are shown, and only every second synthesized beam. 



subtracted down to 5 mJy/beam. For Q, U and /, data 
from all five frequencies are shown in the same plot. The 
RM plot of the Auriga region contains the distribution of 
observed RM (dashed line), as well as that of RM where 
the best-fit linear gradient in RM has been subtracted 
(solid line). The statistical information along these two 
separate lines of sight will be used in Sect.|H|to infer in- 
formation on the Galactic magnetic field and correlation 
lengths in the warm ionized ISM. 

3. Depth depolarization 

The absence of correlated structure in / appears to be a 
general feature: both regions discussed here show it, as 
do other observations made with the WSRT at frequen- 
cies around 350 MHz (e.g. Katgert & de Bruyn 1999, 
Schnitzeler et al., in prep). The lack of corresponding 
structure in / suggests that Faraday rotation is the main 
process responsible for the observed structure in polariza- 
tion. 

However, rotation of the polarization angle cannot by 
itself cause structure in P. As shown in Haverkorn et al. 
(2004), the structure in P in the observations discussed 
here cannot be caused by a missing large-scale structure 
component. Instead, depolarization must be the main cre- 
ator of fluctuations in P. Depolarization can essentially 



occur in three ways: along the line of sight, in the plane 
of the sky or within the frequency band. As the latter 
process only causes significant depolarization for band- 
widths much wider than those used here, we will ignore 
bandwidth depolarization. Depolarization along the line 
of sight is referred to here as depth depolarization (a 
combination of internal Faraday dispersion and differen- 
tial Faraday rotation (Sokoloff et al. 1998, Fletcher et al. 
2004)), as it occurs due to averaging of polarization an- 
gles along the line of sight. In addition, a finite width of 
the telescope beam can cause beam depolarization if the 
structure in polarization angle is on scales smaller than 
the beam. Beam depolarization is clearly observed in the 
fields discussed here, but it cannot explain the structure 
on scales larger than that of the beam (Haverkorn et al. 
2004). 

We are therefore led to consider the situation in which 
the medium that produces the Faraday rotation is mixed 
with a medium that emits synchrotron radiation, which 
produces structure in P through depth depolarization. 
The low level of small-scale structure in / in all of these 
observations implies that the total intensity of the syn- 
chrotron emission must be very uniform. However, the 
magnetic field, and therefore the synchrotron emissivity 
of the thin disk, is far from uniform (e.g. Beck et al. 1996, 
Han & Wielebinski 2002). Therefore, the smoothness of 
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the synchrotron total intensity cannot be due to homoge- 
neous synchrotron emission. Instead, the number of tur- 
bulent cells along the line of sight must be so large that 
the spatial variation in synchrotron emissivity, which is 
a scalar quantity, is averaged out. Linear polarization is 
a vector, so that small-scale structure is more easily pre- 
served. Furthermore, total intensity is integrated over a 
much larger path length than polarized intensity because 
it isn't depolarized. Finally, in the second quadrant, where 
our observations were done, the perpendicular component 
of the uniform Galactic magnetic field is believed to domi- 
nate the component parallel to the line of sight (e.g. Beck 
2001). This means that B±, and therefore the emitted 
synchrotron radiation, has a large uniform component. 

4. Relevant components of the ISM 

4.1. Cosmic rays and thermal gas 

The synchrotron intensity is the integrated non-thermal 
emission along the line of sight. The intensity depends on 
the relativistic electron density Ue^rei and magnetic field 
perpendicular to the line of sight _Bx- Beuermann et al. 
(1985) have modeled the Galactic synchrotron emissivity 
e from the continuum survey by Haslam et al. (1981, 1982) 
at 408 MHz. They incorporate spiral structure in the syn- 
chrotron radiation and find two components of emission, 
viz. a galactocentric thick and thin disk with half equiva- 
lent widths of h^^i, Ri 1800 pc and ft-e,„ « 180 pc at the ra- 
dius of the Sun, respectively, for a galactocentric radius of 
the Sun Rq = 10 kpc. This corresponds to an exponential 
scale height of hsyn,thick ~ 1500 pc and hsyn.thm ~ 150 pc 
scaled to a galactocentric radius of the Sun i?© — 8.5 kpc 
(Beck 2001). Note that the scale heights of the cosmic-ray 
electrons and of the magnetic field must be larger than 
that of the synchrotron disk, e.g. by factors 2 and 4 in 
case of equipartition (Beck 2001). 

The major part of the Faraday rotation is caused 
by the warm ionized medium, contained in the so-called 
Reynolds layer (Reynolds 1989, 1991). The layer has a 
height of about 1 kpc, a temperature T sa 8000 K, 
and a thermal electron density concentrated in clumps of 
ne,th ~ 0.08 cm^'^ with a filling factor of about 20%. 

In this paper we will consider two domains in the ISM. 
The first domain is the thin synchrotron disk, which coin- 
cides with the stellar disk {^-^ 200 — 300 pc), the thin Hi 
disk (~ 200 pc, Dickey & Lockman 1990), and the disk of 
classical H ii regions {^ 60 pc). The second domain is the 
Reynolds layer, which coincides with the thick synchrotron 
disk. There is depolarization in both domains. The Local 
Bubble and Local Interstellar cloud contributions to the 
RM are so small that they are neglected here. 

4.2. Regular and random Galactic magnetic field 

We decompose the Galactic magnetic field in a reg- 
ular large-scale component and a random component 
B = Hreg + '^ran- Estimates of the ratio of random to reg- 



ular magnetic field strengths Bran/ B^eg in the literature 
seem to depend on the method used. Magnetic field de- 
terminations using RMs from extragalactic sources yield 
Bran/Breg ~ 0.5—1 (Jokipii & Lcrchc 1969, Clegg et 
al. 1992). Pulsar RMs indicate that Bran/Breg « 3 - 4 
(Rand & Kulkarni 1989, Ohno & Shibata 1993), although 
this value may be an overestimate (Heiles 1996, Beck et al. 
2003). From diffuse polarization measurements, Spoelstra 



(1984) estimates Bran/Br 



1 — 3, in agreement with 



PhiUipps et al. (1981) who find that Bran/Breg ^ 1- 
Heiles (1996) estimates an average from several studies 

as ^ranj ^reg ^^ i.D. 

Structure in RM is estimated to be present on scales at 
least from 0.1 to 100 pc from observations of extragalac- 
tic point sources (Clegg et al. 1992, Minter & Spangler 
1996), whereas pulsar RMs and dispersion measures DM 
give cell sizes of 10 to 100 pc (Ohno & Shibata 1993). Beck 
et al. (1999) found scale sizes of ^^20 pc for the galaxy 
NGC 6946. 

Field strengths and structure in the Galactic halo, i.e. 
in the gas above the thin synchrotron disk at /i ^ 200 pc, 
can be estimated from observations of halos of external 
galaxies. In observations of synchrotron emission in ha- 
los of edge-on galaxies, the degree of polarization mostly 
increases with distance from the plane, suggesting a de- 
creasing irregular magnetic field component for increasing 
distance to the plane of the galaxy. Structure in the halo 
varies on much larger scales than in the thin disk, viz. on 
scales of about 100 - 1000 pc (e.g. Dumke et al. 1995). 

5. A model of a Faraday-rotating and 
synchrotron-emitting layer 

In this section we describe a simple model of a thin 
Galactic disk containing cosmic rays, magnetic fields and 
thermal electrons, irradiated by a uniform polarized back- 
ground from the thick synchrotron disk. We calculate the 
total intensity, Stokes Q and [/, and the implied RM , for 
various assumptions about the structure of the layer. In 
Sect. El we will compare these results with the observa- 
tions. 

5.1. Outline of the model 

We describe structure in the warm gas and in the mag- 
netic field in the thin disk with a single-cell-size model. 
Fig. 121 gives a sketch of the model and its parameters; in 
addition to the cell size d these are _D, the vertical thick- 
ness of the layer, and the synchrotron emissivity Ic in each 
cell. The warm ionized medium has a filling factor /; this 
is accounted for in the model by setting the electron den- 
sity to an assumed constant value Ue in a fraction / of the 
cells along the line of sight, which are randomly chosen. In 
the remaining fraction (1 — /) of cells, Ue is set to zero as 
an approximation for both the hot dilute gas and the cold 
neutral medium. Thus, we have made the simplifying as- 
sumption that neither the hot nor the cold gas contribute 
significantly to the RM . 
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Fig. 3. Depth depolarization model on a grid of thickness 
Z), containing cells with cell size d. Synchrotron emission 
/c is emitted in each cell, Faraday rotation only occurs in a 
fraction / of the cells. A constant background polarization 
Pfo is also Faraday-rotated while propagating through the 
medium. Along 2 lines of sight (Auriga and Horologium) 
the model is compared to observations. 



with < 77fc < 1 . Details of the emission and propagation 
of the radiation can be found in Appendix fXI 

Each line of sight through the model grid (to be iden- 
tified with the direction of one of our fields, and corre- 
sponding to a particular Galactic latitude) is simulated 
many times, by independently filling the cells that con- 
tain the warm ISM, and by redrawing the angle that the 
random component of the magnetic field makes with the 
line of sight. An ensemble of such realizations, for which 
we derive the distributions of /, Q, U and RM^ simulates 
the many lines of sight for which we obtain the same in- 
formation in one of the observed fields. So, only the statis- 
tical information of these distributions is included in the 
model, not the discrete structure. Each of the two observed 
regions separately already provides useful constraints for 
the model parameters (see below), but the combined set 
of constraints in Table |2] is quite powerful by virtue of the 
different path lengths through the medium. 

As we model many lines of sight by redrawing the same 
line of sight many times for different models, beam depo- 
larization is not included in the models. 



The magnetic field in the thin turbulent disk consists ^'^^ ^^^ ^^^"""^ ^yP^^ °^ '"°^^' parameters 



of a random and a regular component Bj-an and B^eg ■ The 
field strengths of both components are assumed constant 
(but not equal) . As it is not known how B^an and B^eg in 
the cold, warm and hot phases of the ISM are related, we 
consider two extreme cases: 



A: properties of both the random and the regular compo- 
nent of the magnetic field are identical in all phases of 
the ISM. 

B: the random component of the magnetic field only ex- 
ists in the turbulent warm ISM. In the cold and hot 
ISM, the regular magnetic field component predomi- 
nates. The total magnetic field energy density is equal 
in all phases. 

Hence, the properties of the cells are identical in both 
models - they contain ng, B^an and B^eg- The difference 
between models A and B is the hot/cold medium (outside 
the cells) which contains B^an and B^^g (model A), but 
only Breg in model B. 

In each cell, an amount of synchrotron radiation I^ oc 
B\ is emitted, which is assumed to be 70% linearly polar- 
ized. In the cell, the polarized component of this emission 
and all emission from behind is Faraday rotated by an 
amount (^ft- 

The thick synchrotron disk serves as a background to 
the detailed model of the thin disk described here. We as- 
sume that the structure in the thick disk is on such large 
scales that we can approximate the background as a uni- 
form synchrotron emitter, producing a constant polarized 
intensity P\y as input to the thin disk, with uniform polar- 
ization angle. P\, = 0.7 rih lb, where /;, is the total inten- 
sity of the background, and rji, is a constant depolarization 
factor describing the depolarization due to the thick disk. 



Four types of parameters, listed in Table r^TTl are used in 
the model. We discuss these separately. 

Input parameters with fixed values These are physical pa- 
rameters which can be estimated or for which reason- 
ably good estimates exist in the literature. From disper- 
sion measures (DM) of pulsars in globular clusters at 
high Galactic latitude and Ha emission measures (EM), 
Reynolds (1991) derives rig ~ 0.08 cm"-^, with a filling fac- 
tor / = 40% if the warm ionized ISM layer has a constant 
electron density, and 20% if the electron density distribu- 
tion is exponential. The Bcucrmann ct al. (1985) model for 
Galactic synchrotron radiation predicts a half equivalent 
width of the thin disk of 180 pc. We run the model with 
fixed values / = 20%, D = 180 pc and Ue = 0.08 cm"^ 
and discuss afterwards how the results would change if 
these parameters were different. 

The intrinsic polarization angle of the background only 
defines the average angle in the final map of polarization 
angle. It changes the Q and U maps locally, but has no 
influence on the distributions of Q and U. Therefore the 
value of 00 is arbitrary and chosen to be 0°. The angles 
a and 0^ define the orientation of the random component 
of the magnetic field, with respect to the line of sight and 
some fixed direction in the plane of the sky, respectively. 
Both are randomly drawn from uniform distributions, for 
each cell. 

The total intensity Iq is taken from the 408 MHz all- 
sky survey by Haslam et al. (1982). The 2.7 K contribution 
from the CMBR is subtracted before these data are con- 
verted to a frequency of 350 MHz using a spectral index of 
-2.7. Approximately 25% of the total background temper- 
ature is due to point sources (from source counts. Bridle 
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Parameter 


Auriga region 


Horologium region 


RMo 


-3.4 rad m 


-2 


-1.4 rad m~^ 


CTrm 


1.8 rad m" 


-2 


4.3 rad m'^ 


o-/ 


< 1.7 K 




< 2.5 K 


'^Q,t/ 


3K 




4.3 K 


Additional constraints: 






Distributions 


of RM, /, Q, 


and U are Gaussian 



Parameter depends on 



Table 2. Values of observationally determined parameters 
from polarization maps of the Auriga and Horologium re- 
gions and other observational constraints for the models. 



et al. 1972), so only the remaining 75% is included in the 
model. This yields values of 34 K and 47 K for /q in Auriga 
and Horologium, respectively. 

The proportionality constant C = 1 is estimated from 
the local cosmic ray spectrum, see Appendix ^ If strict 
equipartition between cosmic rays and magnetic field ap- 
plies, then C is not constant but varies with B\, so that 
the synchrotron emission / ex B\. Although equiparti- 
tion is believed to hold on global Galactic scales, it is 
highly uncertain if equipartition is valid at parsec scales 
as well, because fluctuations in the supply rate of cosmic 
rays may destroy equipartition on small scales. Therefore, 
the exponent a of the relation / ex i?" could be between 
2 and 4. Although Burn (1966) concludes that equiparti- 
tion does not influence the depolarization much, Sokoloff 
et al. (1998) find that in the case of equipartition the depo- 
larization effects can differ by maximally 25%. We assume 
a := 2 in the model, and discuss afterwards the change in 
parameter values if a > 2. 

Free input parameters No external constraints are imposed 
on the cell size d. Estimates of cell sizes from the literature 
range from 10 pc to about 100 pc, and mostly probe scales 
that exceed the size of our fields. We probe cell sizes from 
a parsec to several tens of parsecs, and find the cell size 
determined in a reasonably narrow range because of the 
observational constraints. 



Constraints determined from the observations As discussed 
in Sect. [3 our observations yield distributions of /, Q, 
U and RM . We summarize these results in Table |21 in 
the form of the mean value of observed RM , RMq, and 
the dispersions UflM, o"/, ag and ajj. In the Auriga field, 
the dispersion in RM was derived after subtraction of the 
best-fit gradient in RM (see Sect.EJ- Like in the observa- 
tions, only reliably determined RMs are used. 

Model parameters that can be adjusted and optimized For 
any chosen value of cell size d, the model contains five pa- 
rameters that are not derived from external data or from 
the observations. These are: the parallel and perpendic- 
ular components of the large-scale magnetic field, Bj.eg,\\ 
and Breg,± respectively, the strength of the random com- 



RMo 


Breg,\ 


1 




C^RM 


Breg,\ 


1 , Bran 




o-/ 


J^ran: 


Breg,± 




<^Q,U 


J^ran^ 


Breg,±, 


Pb 


h 


^ran^ 


Breg,±, 


Pb, rit 



Table 3. Observables and their dependencies 



ponent of the magnetic field Bran, the intensity of the po- 
larized emission from the thick disk Pb, and the thick disk 
depolarization factor rjt, which connects Pf, and /&. The five 
parameters have specific dependences on the observables, 
as is depicted in Tabled This makes it possible to deter- 
mine definite values for most free parameters, e.g. Preg,|| 
can be determined because RMq only depends on Breg,\\- 
The free parameters are determined as follows (where the 
subscript ohs means the observed value): 



Set PregJi to obtain PA/q = PMq, 



ohs 



2. Set Bran to obtain CTrm = o-RM^obs 

3. Set Breg±^ to match ai ~ (Ti,obs- Then: 

- If (jQ^u < crQ,u,obs'- set Pb > to obtain aq^u = 

(^Q,UMbs 

— li<jQu > <jQ^u,obs- decrease Breg^ and find a range 
of {Breg^, Pb) for which ai < aj^obs and aq^u = 
(jQUobs, with the additional constraint that the Q 
and U distributions remain Gaussian. 

4. Set r^b to obtain the correct value of /q 

An example of the correlations used is given in Fig. 01 
The leftmost plot shows the dependence of RMq on P^eg, || • 
As expected, a large regular magnetic field component 
causes a large non-zero mean RM (the negative value of 
the magnetic field reflects the observed negative RMq). 
The center plot gives the change of CTrm with Bran for 
three fixed values of Breg,\\, which increases roughly lin- 
early with Bran and shows hardly any dependence of CTrm 
on Breg,\\- Having set Bran to obtain the observed value 
of CTrm the right plot shows how the observable tqjj de- 
pends on Breg,±- The width of the Q and U distribution 
depends slightly on the chosen values of Bran- 

6. Results from the model 

For models A and B as defined in Sect. 15.11 the propa- 
gation of polarized radiation through the medium is com- 
puted for a range of values of the cell size d, which results 
in values for the five adjustable parameters for each d. 

The allowed cell size is well-constrained by the obser- 
vations: if the cell size is large, the number of cells is small 
for a given path length and filling factor. As the num- 
ber of cells with Faraday-rotating, thermal medium can 
differ per line of sight, the RM distribution will not be 
Gaussian anymore if the cells are chosen too large. On 
the other hand, if the cell size is small, the RM per cell 
decreases. But to obtain a large enough CTrm, the RM 
per cell has to be rather high, so the parameter Bran has 
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Fig. 4. Dependences of depth depolarization parameters on observables in the Auriga region, for model A. Left plot: 
RMq increases with Breg,\\- Center plot: for fixed Breg.\\, ctrm increases with Bram with only a weak dependence on 
the value of Breg.\\- Right plot: ctq^u depends on Bran and Breg.±- 



to be increased to produce the observed value of CTrm- 
However, an increase of Bran increases u/, which then 
puts an upper limit on Breg,±- To produce the observed 
dispersion in Q and U, we then need a large value for 
the background polarized intensity Pi,. If the cell size is 
taken too small, Pb becomes so large compared to the po- 
larized emission in the cells, that the distributions of Q 
and U become distinctly non-Gaussian. Allowed values of 
d range from approximately 1 to 60 pc, with an optimum 
value of about 15 pc, in good agreement with estimates by 
Ohno & Shibata (1993). However, a cell size of 15 pc lo- 
cated at the far end of the thin disk in the direction of the 
Auriga region subtends an angle of more than a degree on 
the sky. As we observe structure on degree to arcminute 
scales, smaller cells must be present as well. Most likely, 
a power law spectrum of turbulence is present on these 
scales (Clegg et al. 1992, Armstrong et al. 1995, Minter 
&: Spangler 1996). An extension of the model including a 
power law spectrum of cell sizes would increase the depo- 
larization along the line of sight, thereby decreasing the 
random magnetic field strength needed. 

We now summarize the result of the comparison be- 
tween the models and the observations. The cell sizes 
probed in the modeling were 1, 2, 5, 10, 20, 30, 40 and 
60 pc, although for model B, only cell sizes above 5 pc are 
allowed, and smaller cell sizes are allowed for the Auriga 
region than for Horologium. 

Fig. |S1 shows the allowed ranges of parameters in the 
two regions for models A (left) and B (right). The up- 
per plots show values of the obtained parallel regular field 
Breg.\\, where the solid line denotes Auriga and the dot- 
ted line Horologium. Breg,\\ is —0.42 ±0.02 /iG for Auriga 
and —0.085 ± 0.005 /iG for Horologium in model A, and 
-0.35 ± 0.01 /iG and -0.065 ± 0.005 /iG respectively in 
model B, where the errors are estimated from the spread 
in observed values. Breg.n hardly depends on cell size. The 
next plots down show Bram where the best value is about 
1 /iG for large (^ 5 pc) cell sizes for both the Auriga and 
the Horologium region, and increases to ^ 4 /iG for cell 
sizes of a parsec. 

For the remaining parameters Breg,±, Pb, and rjb only 
parameter ranges could be determined instead of definite 
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Fig. 5. Allowed ranges of parameters for model A (left) 
and B (right) for different values of cell sizes. The mag- 
netic field values are given in /iG, the polarized brightness 
temperature of the background Pb in K and lo^thin is the 
percentage of the total emissivity that originates in the 
thin disk. The lines in the upper plots show values found 
for Breg.w and Bran for the Auriga region (solid hue) and 
the Horologium region (dotted line). All plots below those 
show allowed ranges in parameters for the Auriga region 
(dark) and Horologium (light). 



values, given in Fig. by a dark grey region for Auriga 
and light grey for Horologium. The edges of the shaded 
areas are not sharp as drawn but have a considerable un- 
certainty. The allowed parameter range should be read 
more as an indication of possible parameter values than 
as ranges with sharp boundaries. Moreover, the discrete 
edges at a certain cell size only indicate that the next 
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probed cell size could not produce parameters in agree- 
ment with the observations. 

The perpendicular component of the regular magnetic 
field Breg,± IS approximately 2.8 ± 0.5 /iG in Auriga and 
3.2 ± 0.5 /iG in Horologium for both models. The inten- 
sity of the polarized background P(, varies between 0.1 and 
3 K, with a best estimate of about 1.5 ± 1.0 K. The depo- 
larization factor 77ft — Ph/(0.7/(,) of the thick disk ranges 
from almost zero to 0.6 with a best estimate of about 
0.15 ± 0.1. 

Below that, Bran/ B^eg IS glvcn, which varies between 
0.5 and 1.5 but is mostly smaller than one. The bottom 
plots show Io,thin, the percentage of the total synchrotron 
emission generated in the thin disk, to be between 20% 
and 75%, and decreasing for larger cell sizes. 

Having determined these parameter ranges, we vary 
the filling factor /, thermal electron density rig or thick- 
ness D while keeping all earlier determined parameters 
fixed, to gauge their influence on the model output pa- 
rameters. A filling factor / ^ 5 - 10% is not allowed in 
either model: large cell sizes give a non-Gaussian RM dis- 
tribution, and small cell sizes yield too high a background 
polarization to keep Q and U Gaussian. No upper limit 
can be given for the filling factor, and Bj-anl B^eg decreases 
with a factor two for / = 1. For varying thermal electron 
density, a low rig ^ 0.03 cm^'^ dictates such a high Bj-an 
that the ratio Breg,±/Pb becomes so low that Q and U be- 
come distinctly non-Gaussian. High electron densities are 
allowed in the models but the random magnetic field drops 
to very low values {Bran ^ 0.15 /iG for Ue ^ 0.1 cm^^^). 
A lower limit to the thickness of the thin disk is about 
100 pc, again no upper limit can be set. The difference 
between a thickness of 150 pc or 180 pc is negligible. 

We checked the influence of the assumption of equipar- 
tition between magnetic fleld and cosmic rays. If / ex -B^ 
instead of / oc B^ , the upper limit to structure in / be- 
comes much more stringent. Therefore, model A will no 
longer produce any solutions that agree with the observ- 
ables. In model B solutions are found with cell sizes 10 
and 20 pc, and Breg,± much lower, about 0.5 fiG. Other 
parameters are comparable to the case where / ex B^ . 

Finally, it should be mentioned that the mean values of 
the distributions of Q and U , which are large-scale compo- 
nents that are not observable with an interferometer, are 
lower than 1.2 K in all models for all parameters. These 
large-scale components are negligibly small, in agreement 
with Haverkorn et al. (2004). 



6.1. Discussion of the resulting model parameters 

A basic flrst conclusion is that the values obtained in the 
two regions roughly agree, even though the Auriga and 
Horologium regions have different input parameters and 
a different line of sight through the medium. The regular 
magnetic field components approximately agree in the two 
regions {Breg ~ 3 /iG). 



From the deduced depolarization factor rji, = 0.15, we 
can estimate the halo magnetic field, assuming a constant 
background. Around 350 MHz, depth depolarization of a 
uniform medium to 0.15 times the original polarization is 
caused by RAI « 3 — 5 rad m^^ (Burn 1966), which indi- 
cates a value of i3|| w 0.1/iG for a height of the Reynolds 
layer of 1 kpc and rie ~ 0.05 cm~^ in the halo. So the disk 
magnetic field could persist with only slight attenuation 
throughout the Reynolds layer, as was suggested earlier 
by Han et al. (1999). 

Our estimate of Bran/ Breg < 1 is somewhat lower 
than most of the estimates from the literature discussed 
in Section 01 This may be due to several factors. First, 
random magnetic field structure on scales larger than our 
field of view (~ 7° x 9°) will be interpreted as regular field 
in our analysis. Secondly, it could be the result of selection, 
as our observational fields were chosen for their high po- 
larized intensity, which in our model automatically implies 
a modest random magnetic fleld. Finally, the observations 
are in the second Galactic quadrant, so we probe mostly 
the inter-arm region between the Local and the Perseus 
arms, where Bran/ Breg is smaller than in the average ISM 
including spiral arms (e.g. Indrani and Deshpande 1998, 
Beck 2001). 

The emission in the thin disk lo^thin is also estimated 
by Beuermann et al. (1985) in their standard decomposi- 
tion of /q into thin and thick disk contributions. According 
to their model, only about 20 to at most 35% of Jq is 
generated in the thin disk and the nearest 180 pc of 
the thick disk. Furthermore, Caswell (1976) estimated the 
synchrotron emissivity from a survey with the Penticton 
10 MHz array as 240 K pc^^ at 10 MHz. Rescaled to 
350 MHz, this gives a total emission from the thin disk of 
10.6 K in the Auriga region and 21 K in the Horologium 
region. Roger et al. (1999) estimate from the 22 MHz sur- 
vey performed with the DRAO 22 MHz radio telescope 
an emissivity of about 55 K pc~^ for two HII regions in 
the outer Galaxy, out of the Galactic plane. Their results 
give estimates of the emissivity in the thin disk which 
are approximately twice as high as the estimates from the 
Caswell survey. Due to the large uncertainty in the emis- 
sivity, Io,thin does not put a strong constraint on the model 
parameters. 



6.2. A "polarization horizon"? 

A "polarization horizon" is defined as a distance beyond 
which (most of the) emitted polarized emission is depolar- 
ized when it reaches the observer (Landecker et al. 2001). 
This can be due to beam depolarization, when the angular 
scale of the structure in the polarized emission becomes 
smaller than the synthesized beam at a certain distance. If 
the smallest scales in the observed regions are about a par- 
sec, the angle of these scales on the sky becomes smaller 
than the beam at a distance of about 700 pc. Spoelstra 
(1984) derived the polarization horizon from comparison 
of radio continuum data at five frequencies from 408 MHz 
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Fig. 6. Observed fraction of polarization originating closer 
than a certain distance against the path length, for the 
Auriga region (solid line) and the Horologium region 
(dashed line). 



to 1411 MHz (Brouw & Spoelstra 1976) with starhght 
polarization. He found a distance to the origin of the po- 
larized radio emission of 625 ± 125 pc in the direction of 
our fields. 

Furthermore, the resulting degree of polarization de- 
creases for increasing pathlength through a rotating and 
emitting medium. However, as depolarization is a process 
occurring in a telescope and not in the medium itself, it 
cannot be determined what distance the remaining polar- 
ized radiation originated. We can only estimate the de- 
crease in polarization as a function of path length. In our 
model, we build up a line of sight by adding cells one 
by one, starting at the observer. The radiation from each 
added cell is Faraday-rotated by all warm ionized mate- 
rial in front of it. The observed degree of polarization after 
addition of each cell is given as a function of the line of 
sight built up until that particular cell in Fig. El for the 
Auriga (solid line) and Horologium (dashed line) regions, 
for model A with a cell size of 5 pc. Even for the to- 
tal path length through the thin disk, still a fairly large 
fraction (20%) of the polarized emission can be observed. 
Therefore, depth depolarization alone cannot produce a 
true horizon, but attenuates the polarization gradually 
with path length. 

From both arguments, we estimate a distance of about 
600 to 700 pc as the critical path length ("polarization 
horizon" ) . Polarized radiation traveling along significantly 
larger path lengths than this is expected to be largely 
depolarized. 



We modeled the effect of depth depolarization with a 
simple model of the Galactic ISM consisting of a layer of 
cells containing random and regular magnetic field Bran 
and Breg, and thermal electron density Ue in a fraction / 
of the cells (mimicking the filling factor /). This layer cor- 
responds to the Galactic thin disk with small-scale struc- 
ture in the magnetic field. The Galactic thick disk or halo 
is modeled by a constant background Pb, with a certain 
constant depolarization denoted by the factor rjb. We vary 
cell size, magnetic field and background to obtain a range 
of models that comply to the observational constraints, 
i.e. yield the correct width, mean and shape of the distri- 
butions of Q, U , I and RM. 

The results can be summarized as follows: the allowed 
cell size is constrained to be in the range of 1 to 60 pc, 
with a best estimate of 15 pc. The regular magnetic field 
component along the line of sight (~ —0.4 /iG for Auriga, 
and ^ —0.08 fiG for Horologium) is much smaller than 
the regular magnetic field component perpendicular to the 
line of sight (^ 2.8 /iG for Auriga, and ^ 3.2 /iG for 
Horologium), indicating that the regular magnetic field is 
directed almost perpendicular to the line of sight in these 
directions. The random magnetic field component is about 
1 to 3 /iG in the two regions. In most of our models, the 
regular component of the magnetic field was found to be 
higher than the random component, with an average ra- 
tio of Bran/ Breg — 0.7 ± 0.5 which incrcascs for smaller 
cell sizes. Estimates from the literature tend to give larger 
ratios (0.5 ~ 4). This could be explained by the size of 
our fields of view (~ 7° in size) , so that random compo- 
nents of the magnetic field on scales large than the field 
are misinterpreted as regular components. Furthermore, 
the fields of observation were selected for their high polar- 
ization, indicating a higher regular magnetic field compo- 
nent than average, and are situated in an inter-arm region, 
where uniform magnetic fields tend to be higher than av- 
erage. The constant polarized background intensity from 
the thick disk is about 1.5 ± 1.0 K. 

This model forms a promising first attempt to derive 
properties of the Galactic magnetic field from observed 
polarization and rotation measures. In future work, the 
model can be expanded e.g. using a power law distribu- 
tion of the structure. Furthermore, cell size, filling factor 
and electron density appear to be correlated (Berkhuijsen 
1999), which should be incorporated in a future version. 
New observations in different directions can narrow down 
the parameter space considerably. 



7. Summary and conclusions 

Depth depolarization, the depolarization process along the 
line of sight in a medium of synchrotron-emitting and 
Faraday-rotating material, is the dominant cause of struc- 
ture in polarized intensity which is unrelated to total in- 
tensity fluctuations. 
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Appendix A: Outline of the depth depolarization 
model 

The synchrotron radiation emitted in each cell is Ic oc B\ . 
This emission, and the emission from each cell further 
away along the line of sight and from the background pass- 
ing through the cell, is Faraday rotated by an amount cjipr- 
So in each cell: 



Ic 
Pr 



[{Branii'^^af + Bl^gj] 



C 

— [(^i(™„smaj'-+±(;^g^ 

RM\^ = 0.81 Tie {Bran COS a + Breg,\\) d\^ 
0.7 Ir 



(A.l) 



where Bran is the constant strength of the random mag- 
netic field in /iG, a its random angle with the line of sight, 
C a proportionality constant, N the number of cells along 
the line of sight and d the path length. The total emis- 
sion from the layer (« CN) is comparable for different 
cell sizes, therefore a factor 1/A^ is added to Eq. (|A.1|I . 
The polarized emission in each cell Pc equals the maxi- 
mum polarization of synchrotron radiation Ic generated 
in a cell. For an electron energy power law distribution 
N{E) ex E ', the degree of polarization p is related to 
the spectral index 7 of the electron energy distribution as 
p(7) = (37 + 3)/(37 + 7) (Burn 1966). For 7 around 2.7, 
the maximum polarization is ^ 70% of the total intensity. 
The polarization angle of the emission generated in each 
cell (j)in is taken to be perpendicular to the position angle 
of the perpendicular magnetic field. The position angle of 
the random magnetic field component (f>r is random, and 
that for the regular component is chosen in the direction 
of Galactic longitude. Therefore the polarized intensity 
emerging fronr a cell is 



Pc = 0.7/e e 



-2i{4>Fr+cl>z^) 



0.7/6 e" 



-2l(<t,Fr-+lt>b) 



(A.2) 



for a cell that is irradiated with polarized intensity lb and 
polarization angle (j)b. The input and output parameters 
are given in Table IXtI 
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Input parameters with fixed values 


value: 


He thermal electron density in cells 


0.08 cm--^' (Reynolds 1991) 




/ filling factor of the warm ISM 


20% (Reynolds 1991) 




D thickness of the layer with cells 


180 pc (Beuermann et al. 1985) 




(f>o intrinsic polarization angle of the background 


arbitrary: 0° chosen 




(f>r position angle of random magnetic field 


random per cell 




a angle between random magnetic field and line of sight 


random per cell 




Jo total intensity 


Auriga: 34 K 






Horologium: 47 K (from Haslam et al. 


1982) 


C proportionality constant between 7c and B^ 


C = 1, see Appendix ini 




Free input parameters 


d cell size 


Constraints determined from the observations 


RMo mean rotation measure 






ctrm width of RM distribution 






ai width of / distribution 






<^Q.u width oi Q, U distribution 






Model parameters that can be adjusted and optimized 


set by dependence of: | 


-Brea.ii parallel component of regular magnetic field 


RMo (Bres.ii) 




Bran (constaut) strength of random magnetic field 


0"rm (-Breg, II , -Bran) 




Breg.± perpendicular component of regular magnetic field 


O"/ (Bran, Breg, ±) 




Pi polarized intensity of background 


<^Q,U (Bran, Breg, ±,Pb) 




rjb factor for depolarization of background 


/() [Bran, Breg, ±,Pb,'nb) 




Additional constraints: 


background depolarization factor < 775 < 1 






Number of cells A'^ — L/d, while Nf cells determine the shape of RM distribution 




Breg,±/Pb determines shape of Q, U distribution 







Table A.l. The first set of parameters is determined from the literature or can be arbitrarily chosen. The second set 
is varied in the models, and the third set of parameters is set by our observations. The last set are those parameters 
of the ISM that can be estimated from the models, followed by the input parameters from the categories above. In 
parentheses the model parameters that they depend on. 



Appendix B: Estimate of the parameter C 

The total power per unit volume per unit frequency of 
synchrotron emission is (Rybicki & Lightman 1979) 



Ptotiio) 



Viq^n'B sm( 
27r?Tic^(p + 1 

/ mcuj 

X : 

\ 3qB sin a 

3V3q^K'Bl_ 



rf^ + i^ 

V4 12 

-(p-l)/2 



4 12 



function of energy N{E)dE = KE~PdE derived from di- 
rect measurements of the particle spectrum in the local 
ISM, in agreement with the value found by Golden et al. 
(1994). This can be converted into k' as k! — K{rac^)'^^~P\ 
Using p — Z and converting to cgs units gives 
1.0 10^-* 
yields 



part cm ^ ergs^. Inserting this into Eq. 



forp = 3 (B.l) ^totM=4.2 10 



-39 



J/ [MHz] 



W m--^ Hz- 



(B.2) 



where q is the electron charge, k' is the proportionality 
constant in a power law particle spectrum iV(7) with spec- 
tral index p {N{j)d'y — n'j^Pd'y), Bsina = B± is the 
magnetic field component perpendicular to line of sight, 
m is the electron mass, c the speed of light and uj the 
angular frequency of the radiation. T denotes the gamma- 
function. 

We can estimate k' by assuming that the electron par- 
ticle spectrum throughout the ISM is equal to the local 
value in the solar neighborhood. Longair (1981) gives a 
value of K = 2.9 10"^ particles m'^ GeV-^^-^) for the 
proportionality constant of the particle spectrum as a 



This is the volume emissivity of synchrotron emission 
as a function of frequency and magnetic field. We can 
check how reasonable this number is by comparing to 
the total power observed from the Caswell (1976) ra- 
dio survey at 10 MHz. The average brightness temper- 
ature computed by Caswell corresponds to a volume 
emissivity of ~ 3 10"^^ W m'^ Hz"! (Longair 1981), 
which agrees well with our value of Ptot(10 MHz) = 
0.42 10-39 Bj_[nG]^ W m-3 Hz^i for Bj_^3 ^G. 

The next step is to describe Ptot in terms of the ob- 
servables. The observables are in Kelvin, whereas Ptot is 
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in W m"'^ Hz^^. The power that is detected on 1 m^ of 
antenna surface is 

Pceiuant [W Hz-^ m'^] = -^^Ptot [W m-^ Hz-1] 

with D the distance to the cell. This gives the emissivity 
of a source in the direction of the observer on 1 m^ of 
dish in Jansky's. Then, to convert into Jansky/beam, only 
consider the part of the source which fits into a beam. 
With a spatial resolution of Z3tan(5'), the area of a beam 
is 7r(l/2 £'tan(5')^), and the number of beams that fits 
into the source is AcP /{ttD^ tan^(5')). Then the emissivity 
in one cell per beam is 



Ic[W Hz-i m-2 beam-i] = 

Ptot [W m-3 Hz 



d' f7:DHan^i5'),^ rw ^-3 tt^-H 



47rD2 V 4d2 



so that, using Eq. (|B.2|) . the emissivity Ic in Jy/beam is 
Ic = 0.21 d[pc] K using that 1 mJy/beam « 0.13 K at 
350 MHz. Combining this result with Eq. (jA.ip yields an 
estimate for C: 

V 

for a path length L = 900 pc and v = 350 MHz, where we 
used N = L/d. We use a value of C = 1 in the model. 



